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Abstract 

We study the spectrum of rotating Bose-Einstein Condensates in the limit of 
weak repulsive interactions and present analytical and numerical results for 
energies and wave functions. While the low- lying states are of collective na- 
ture, the high-lying states are dominated by single-particle excitations and can 
accurately approximated by simple polynomial expressions. In the limit that 
the number of particles is large compared to the number of excited quanta, 
the single-particle states become excellent approximations to the eigenstates, 
and a rather simple ordering scheme is obtained. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db 



1 



Typeset using REVTgX 



The experimental realization of Bose-Einstein condensation in atomic vapors |T]-g] has 
received much interest in recent years [f|]. Very recently experiments have started to in- 
vestigate Bose-Einstein condensates under rotation. Matthews et al. created vortices in a 
two-component condensate ||, and Madison et al. studied rotations of a one-component 
gas. Theoretical studies have focused on the Thomas-Fermi regime of strong interactions 
f^H or on the limit of weak interactions between the atoms [10-|17|. While the former case 



is closer related to the experiment, the latter case allows for exact (numerical or analytical) 
solutions and is of particular theoretical interest. Wilkin et al. studied the case of attrac- 
tive interactions [TU] and employed a composite boson/fermion picture [12,13] to describe 



configurations beyond the one- vortex state in few-body boson systems. Mottelson [11] and 



the Copenhagen group |l5[|l6j analyzed the mean- field theory in the weak interaction limit. 



In Ref. []T3[ the authors studied the ground states using exact diagonalization techniques. 
So far, research has been restricted to the low-lying states in rotating Bose-Einstein con- 
densates. It is the purpose of this paper to study higher excited states in the rotational 
spectrum. 

The Hamiltonian is 



j=i 1 z z > i<j 



where the trap frequency is set to hu = 1 and g is a dimensionless coupling. In the limit of 
weak repulsive interactions < g <C 1 and maximal alignment of angular momentum with 
the z-axis the Hamiltonian becomes essentially two-dimensional [jn|,[nj . In what follows we 
restrict ourselves to this limit and fix the total angular momentum to L. 

The problem then consists in diagonalizing the contact interaction in the Hilbert space 
of degenerate oscillator states at energy total L (setting the ground state energy to zero) 
p0| . We will analyze the problem in second quantization and in configuration space. While 



analytical results will be derived in both representations, numerical work is most easily 



done using second quantization [18|.[19(]. Let us start with the second quantized two-body 
interaction 

V = y^27r5(rj — r,-) % , f , j=f=j= 5lM a]a\a k ai. (1) 

£rf v JJ 2 . z t, 2 k+l Ji j k l K J 

■i<j v J 

The operators a m and annihilate and create one boson in the single-particle oscillator 
state \m) with energy mhuj and angular momentum mh, respectively and fulfill bosonic 
commutation rules. For total angular momentum L the Fock space is spanned by states |a) = 
|n ,ni, ...,n k ) with J2i=o,k n i = N , a]aj\n ,ni, ...,n k )= nj\n ,ni, ...,n k ) and T,j=o,k3 n j = 
L. Here rij denotes the occupation of the j th single particle state 

Turning to configuration space we use complex single-particle coordinates z = x + iy 
such that (z\j) = <j>j(z) = (7rj!) -1 / 2 z^ exp (—^\z\ 2 ). Basis functions for the many-body 
problem are ip(zi, . . . , z^) rL=i ex P ( — il z j'| 2 ) where ip(zi, . . . , zn) is a homogeneous polyno- 
mial of degree L that is totally symmetric in its coordinates. For representations of such 
polynomials see, e.g. ref. f2lfl . In what follows we will omit the exponentials when working 
in configuration space and identify wave functions directly by their polynomial part. Let 
ip(zi, . . . , Z]\r) be such a polynomial. It is easily seen that the two-body operator 
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2 k+l n\m\k\l\ 



Tk+l 
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z n z m 



dtdl 



(2) 



Z<i — Zj 



has the same matrix elements as the contact interaction 2ir5(ri — r^). Here it is understood 
that the derivatives act only onto the polynomial part of the wave function and that they 
are evaluated at = Zj = 0. Performing the sums in eq. (g) and introducing coordinates 
R%j — { z i + Zj)/ \/2 and = (z f — Zj)/^/2 leads to the following 

27T(5(r i - rj) tp(zx, ...,Zi,...,Zj,...,z N ) = Wij ip{zt, . . . , z u . . . , Zj, . . . , z N ) 

= ip (zi, ..., Zl ^ Zj ,..., Z3 ^ Z \ ...,z N ^j. (3) 

Thus, the contact interaction 2ir5(ri — rj) acts on a wave function by replacing the i th and 
jth ar g um ent by its mean value. Note that the center of mass 

1 N 

is left invariant under the action of the contact interaction. 

Let us now consider the spectrum of states at angular momentum L. States may be 
classified by their angular momentum quantum numbers L and L cm , and by their energies 
E n = L + ge n . The angular momentum of the center of mass and is denoted as L cm and 
given by 



A' 



1 



3=1 7V k I 



(4) 



where the differential operator is understood to act only on the polynomial part of the wave 
function. The functions z r and 



N 



qk(zi, ...,z N )= ^2(zj - z c ) , k > 2 



(5) 



are eigenfunctions of L cm with eigenvalues one and zero, respectively. Products z^q\ 2 . . . q\ N 
with A being a partition of L are eigenstates of L cm with eigenvalue L cm = Ai and span the 
Hilbert space [20| at angular momentum 2 < L < N. Note that they do not form an 
orthogonal basis. Let \L, L cm ,e n ) denote an eigenstate of the Hamiltonian ([I]). Due to an 
5*0(2, 1) symmetry one has z c \L,L cm ,e n ) oc \L + l,L cm + l,e„) p2| , |23[| . It is therefore of 



particular interest to find the "non-spurious" eigenstates (i.e. those with no excitations of 
the center of mass coordinate) since all other states can simply be obtained by multiplication 
with powers of z c |]22|JT0| , [17|| . Figure [l] (left part) shows the spectrum of non-spurious states 
for N = 50 bosons and angular momenta 2 < L < 18. Note that there are several high-lying 
levels that approach constant values as L increases. This is particularly the case for two 
most energetic levels, and these are well separated from their neighbors for all values of L. 
As has been observed earlier in exact diagonalizations |13[ and mean-field calculations |7|,|T5 
the ground states depend linearly on total angular momentum. 
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Based on numerical results, we suggested in ref. |L3| that the ground state wave function 



is 



V>o(*l, ...,Z N )= Yl ( Z Pl - Z c)(z P2 ~ z c) ■ ■ ■ (Z PL - Z c ) 

l<p 1 <p 2 <...<p L <N 



(6) 



and the ground state energy is 



eo = - JV(JV-l)- 



NL 



(7) 



Very recently, Jackson and Kavoulakis [jn[ showed the existence of an eigenstate with energy 
(0) while Smith and Wilkin ]24| proved that the wave function @ indeed has eigenvalue (0) . 
This proof was given using the second quantized form of the interaction. As an application 
of eq. (H) we now present a shorter proof in configuration space. We expand 



vL-fc 



E 

l<pi<...<p k <N 



Zpi • • ■ z Pk . 



Applying operator (fj) to this expansion and using the prescription ([|) yields 

L (i 

k=0 \- 



A fN-k s 
L-k , 



\L-k 



Z pi . . . z Pk 



E 

l< Pl <...<p k <N 

E 

l<pi<...<p fc _i<iV 



+ {Zi + Zj 



+ ^K z i + z j) 



E 



Z pi • • ■ Z Vk-\ 



z Pl ■ ■ ■ z Pk _ 2 



l<Pi<...<p k - 2 <N ) 

Using (zi + Zj) 2 = (zi — Zj) 2 + AziZj in the third term enclosed in brackets one obtains 



Wij ip = ipo + 

k=0 



L-k 



—z, 



\L-k 



^\ Z i Z ji 



E 

l<pi<...<p fe _2<A r 



Z P1 ■ ■ ■ Z Pk-2 



Summing over J2i<j an d using the identity 



E 

Kj 

—kN 



Zi Zj / 



E 

l<pi<...<p fe _2<A r 



Z Pl ■ ■ ■ Z Pk-2 



E z pi ■ ■ ■ z Pk + N ( N -k + 1) z c %i • • • z Pk-i 

l< Pl <...<p k <N l<Pi<-<Pk-i<N 



one arrives directly at the desired result Vipo = eo^o- 

Another exact result is readily obtained for L = 4 and arbitrary N. Employing eq. @ 
one may show that 



v E 



. Zi Zi 



l<i<j<N 



N 
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-N + 



E 



(8) 



l<i<j<N 
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Unfortunately, we are not able to generalize this result to higher values of angular momen- 
tum. Note however that wave functions (§) and comprise the set of non-spurious states 
for L = 4. 

Let us consider the most energetic non-spurious state. Numerical computations show 
that the ansatz (compare with eq. @) 

ipxiz!, . . . , z N ) = q L (zi, . . . , z N ) (9) 

has almost unit overlap with the most energetic non-spurious eigenstate and is approaching 
the asymptotic eigenvalue ei = h(N — 1)(N — 2) exponentially fast with increasing values of 
angular momentum L. Table | presents numerical values obtained for the overlaps with the 
exact numerical eigenf unction, energies {ipi\V\ipi) , and widths = (ipi\V 2 \ipi) — {ipi\V\ipi) 2 . 
The presented results show that eq. (|) is a highly accurate approximation to the exact 
eigenstate. This can be understood in more detail by applying the operator (P to ip\ 

Vik = \{N - 1)(JV - 2) + 2 1 - 1 ]T ( Zi + z 3 - 2z L ) L . (10) 
1 i<j 

For sufficiently large L the "remainder" Yli<j( z i + z j ~ % z c) L is numerically found to have 
significant overlap only with a small number of eigenstates in the bulk. This explains 
the quality of our ansatz. Note that xpo cx ip\ for L = 2 and L = 3. Note also that 
Vipi = [1 + (— l) L 2 1 ~ L }ipi for iV = 3. The ansatz @ describes ipi as an excitation of one 
particle relative to the center of mass. This picture is confirmed by the computation of 
selected single-particle occupation numbers nj = (ipi\nj\ipi) shown in Figure |2|. The main 
contribution to the total angular momentum stems from the single-particle orbital 4>l{ z ) and 
a smaller part is carried by the orbital <Pl_i(z). With view to eq. ([!]) it becomes clear that 
the energy of ip\ depends very weakly on L. The energy contribution of the highly excited 
single-particle states is exponentially suppressed for sufficiently large angular momenta. 

We may also obtain an excellent approximation for the second highest non-spurious state. 
The ansatz 

N — 1 

ip 2 (zi, ...,z N )= — q 2 {z u ...,z N ) qL-2{zi, ■ ■ ■ , z N ) - ^\{z u ...,z N ) (11) 

has almost unit overlap with the numerically obtained exact eigenstate for sufficiently large 
values of total angular momentum. It approaches the numerically determined energy e 2 = 
|[(iV — 2) 2 — 1] exponentially fast with increasing L. Details concerning overlaps and width 
are presented in Table |. The dominant contributions to the total angular momentum of ip 2 
stem from the orbitals <p 2 (z), <Pl-3( z ) and (pL-2( z )- As expected the ansatz (|Tl|) describes 
mainly a few-boson excitation. The single-particle character of the two most energetic states 



differs strongly from the collective behavior found for the the ground state ipo. ||1T| , P^3| . |15 
This observation is similar to the case of bosons in the Thomas-Fermi regime [2^] and to 
the case of many-fermion systems like atomic nuclei p6 |. 



We checked that the numerical results agree with ansatz (|]) and ([11]) for various num- 
bers of bosons. The convergence is determined by L only and depends weakly on N. The 
spectrum depicted in Figure [I] seems to suggest that there are more states which may be 
approximated by simple polynomial expressions; however we do not see an obvious general- 
ization of the expressions (P) and ([TT]) to account for further states. 
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Let us consider the regime iV ^> L ~ 0(N°) where the ground state is occupied by a 
macroscopic number of bosons, i.e. no ~ O(N). While low-lying states have been computed 
in refs. [ IT| , P^ | we aim at a general understanding of the full spectrum within this regime. 



Eq. (0) shows that the diagonal matrix elements of V are of order 0(N 2 ) while off-diagonal 
ones are only of order 0(N 1 ^ 2 ). In a first approximation one may then neglect any off- 
diagonal elements and take the single-particle states as eigenstates. One obtains from eq. ([!]) 

V\n Q ,n h n 2 , . . . ,n k ) = e{ n .y\n ,ni,n 2 , . . . ,n k ) + 0(N 1/2 ), 

e {nj} = -n (n -l) + 2n J2 2 ~ kn k- ( 12 ) 

1 k>0 

Note that typical energies are of the order 0(N 2 ) while energy differences for two different 
states are of order O(N). Small energy differences may occur for two configurations {nj} and 
{n'j} that differ only in the occupation of the orbitals j = 0, 2, 3. Both configurations should 
have identical number of bosons and angular momentum, i.e. hq + n 2 + n% = n' + n 2 + n' 3 
and 2n 2 + 3n3 = 2n 2 + 3n' 3 . The energy difference of such configurations is only of order 
~ 0(N°) which is much smaller than a typical spacing. Thus there are quasi-degenerate 
levels in the spectrum, and the ground state becomes quasi-degenerate for L > 9. This 



observation has also been made by Mottelson |TT 



Note that the single particle states have almost good quantum number L cm . One finds 
L C m\n>o, Tii, n>2, ■ ■ ■ , n>k) ~ ni\n , m, n 2 , . . . , n k ), and corrections are of order 0(L/N). This 
can be seen by applying the operator (f|) to a state with uq ~ 0(N). Thus, the non- 
spurious states are comprised by configurations with n\ = 0. We may now discuss the order 
of the spectrum of non-spurious states. According to eq. ( |T2"D this ordering is led by two 
observations, (i) The lower the number of bosons in the ground state, the lower the total 
energy of this state, (ii) For two configurations with equal numbers of excited bosons the 
higher energetic one has the highest single particle orbital occupied. This ordering scheme 
favors collective motion above single particle excitations. Deviations from this scheme are 
caused by the quasi- degeneracy of configurations discussed above. The validity of our argu- 
ments and of eq. (0) is confirmed by numerically computing and analyzing the spectrum 
for N = 46000 in the regime 2 < L < 18. 

Let us finally examine the iV-dependence of the spectra at fixed L as shown in Figure |3|. 
The spectrum evolves rather smoothly with increasing iV and levels tend to be clustered 
at high N. For instance, there is a triplet of states at lowest energies. This is due to 
the quasi- degeneracy discussed above, and the triplet is comprised by configurations with 
{n 2 = 6}, {n 3 = 4} and {n 2 = 3,n 3 = 2}. Such quasi-degeneracies indicate a high degree of 
integrability and regularity. However, there are also a few quasi-degeneracies at mid-values 
of N. Closer inspection reveals that these are avoided crossings. As an example we present 
in Figure § (top) the avoided crossing at iV w 70. The enlargement of the energy clearly 
indicates an avoided crossing. Further evidence comes from analyzing the wave function 



structure, i.e. we follow the inverse participation [27 



^ = EIW>I 4 (13) 

a 

in the basis of single-particle basis states \a). The inverse 1/i^ measures the number of 
basis states that have significant overlap with Figure |] (bottom) shows that the wave 
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function structure changes at the avoided crossing. The simple ordering obtained for the 
spectrum at large values of N is thus lost when decreasing N. We also recall that avoided 
crossings indicate some degree of level repulsion and chaoticity PB|. Therefore, we do not 



expect that all states may be described by simple polynomial expressions. Note however, 
that the high lying non-spurious states and the ground state do not interact with their 
neighbored levels - a further confirmation of their regularity. 

In summary, we have studied the rotational spectrum of Bose-Einstein condensates in 
the limit of weak repulsive interactions. The most energetic non-spurious states are well 
separated from neighboring levels and their energies approach constant values with increasing 
angular momentum. These states can accurately be approximated by simple polynomial 
expressions. They exhibit single-particle behavior and differ thereby from the collective 
behavior of the ground state. In the limit of a macroscopic occupation of the ground state and 
small total angular momentum the single-particle states become approximate eigenstates. 
We obtain a simple ordering scheme for the entire excitation spectrum that energetically 
favors collective behavior over single-particle behavior. This ordering scheme is expected 
to break down when leaving the limit of L/N <C 1 since some states undergo avoided 
crossings. While this introduces chaotic behavior in the involved states, the ground state 
and the highest excited non-spurious states appear to remain simple and regular. Parts of 
the analytical results presented in this work including a proof concerning the ground state 
wave function are based on a particularly useful and novel representation of the contact 
interaction in configuration space. 
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7 



REFERENCES 



e-mail: papenbro@phys.washington.edu 
e-mail : b er t sch@ phy s . Washington . edu 

M. N. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, 
Science 269, 198 (1995). 

C. C. Bradley, C. A. Sacket, J. J. Toilet, and R. G. Hulet, Phys. Rev. Lett. 75, 1687 
(1995). 

K. B. Davis, M.-O. Mewees, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. 
Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

For a review see, e.g., F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. 
Mod. Phys. 71, 463 (1999). 

M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. Hall, C. E. Wieman, and E. A. 
Cornell, Phys. Rev. Lett. 83 (1999) 2498. 

K. W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, Phys. Rev. Lett. 84, 806 
(2000). 

D. A. Butts and D. S. Rokhsar, Nature 397, 327 (1999). 
M. Linn and A. L. Fetter, Phys. Rev. A 60, 4910 (1999). 

D. L. Feder, C. W. Clark, and B. I. Schneider, Phys. Rev. Lett. 82, 4956 (1999);Phys. 
Rev. A 61, 011601 (2000). 

N. K. Wilkin, J. M. Gunn, and R. A. Smith, Phys. Rev. Lett. 80, 2265 (1998). 
B. Mottelson, Phys. Rev. Lett. 83, 2695 (1999). 
N. R. Cooper and N. K. Wilkin, Phys. Rev. B 60, R16279 (1999). 
G. F. Bertsch and T. Papenbrock, Phys. Rev. Lett. 83, 5412 (1999). 
N. K. Wilkin and J. M. F. Gunn, Phys. Rev. Lett. 84, 6 (2000). 
G. M. Kavoulakis, B. Mottelson, and C. J. Pethick, |arXiv:cond-mat /00043"07| . 



A. D. Jackson, G. M. Kavoulakis, B. Mottelson, and S. M. Reimann, arXivicono 1 "- 
mat/0004309|. 



A. D. Jackson and G. M. Kavoulakis, |arXiv:cond-mat/0005159| . 



T. Haugset and H. Haugerud, Phys. Rev. A 57, 3809 (1998). 

T. Papenbrock and G. F. Bertsch, Phys. Rev. A 58, 4854 (1998). 

G. Natanson and R. S. Berry, Ann. Phys. 155, 158 (1984); Ann. Phys. 155, 178 (1984). 

I. G. Macdonald, Symmetric Functions and Hall Polynomials, 2nd edition, Clarendon 

Press, Oxford 1995. 

L.P Pitaevskii and A. Rosch, Phys. Rev. A 55, R853 (1998). 

A. Perelomov, Generalized Coherent States and Their Applications, (Texts and mono- 
graphs in physics), Springer Verlag, Berlin 1986. 
R. A. Smith and N. K. Wilkin, |arXiv:cond-mat/ 0005230 . 



F. Dalfovo, S. Giorgini, M. Guilleumas, L. Pitaevskii, S. Stringari, Phys. Rev. A 56, 
3840 (1997). 

N. Bohr and B. R. Mottelson, Nuclear Structure, Vol. I (Benjamin, New York, 1969). 
See, e.g. L. Kaplan, Nonlinearity 12, Rl (1999). 

T. Guhr, A. Muller-Groeling, and H. A. Weidenmuller, Phys. Rep. 299, 189 (1998). 



8 



TABLES 



L 


8 


12 


16 


20 


l-KV'llV'exact)! 2 


1176.325 

0.6 
2 x 1(T 4 


1176.019 

0.2 
6 x 10" 6 


1176.001 
0.03 

2 x 10~ 7 


1176.000 
0.01 

< 10- 7 


<Mv\ik) 
r 2 

1- K^IV'exact)! 2 


1153.031 
1.4 

3 x 10~ 3 


1151.604 
0.3 

9 x 10~ 5 


1151.508 
0.09 

5 x 10~ 6 


1151.501 
0.03 

4 x 10~ 7 



TABLE I. Numerical data for the states ipi and ^2 as a function of angular momentum L for 
N = 50. The energy expectation values (^j\V\i(jj) approach energies ei = (N — 1)(JV — 2)/2 = 1176 
and e 2 = ((N - 2) 2 - l)/2 = 1151.5, respectively. The widths L 2 = (^j|L 2 |V>j> - {tpjlV^j) 2 are 



small compared to the level spacing (which is of order O(N)), and the overlap with the exact 
eigenstates is close to one. 
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FIGURES 




FIG. 1. Spectrum of non-spurious states for a system of N = 50 particles as a function of 
angular momentum. 




FIG. 2. Structure of state Single-particle angular momenta j{^i\hj\ipi) as a function of 
angular momentum L for a system of N = 50 bosons, (j = L — 1: circles; j = L: triangles; sum of 
both: diamonds; total angular momentum: dashed line) 
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FIG. 3. Spectrum of non-spurious states for L = 12 as a function on the number of bosons N. 
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FIG. 4. Avoided Crossing. Top: Two energies as a function of N. Bottom: Number of 
participating basis states as a function of energy. Results are taken from one avoided crossing at 
L = 12. 
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